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The NA61/SHINE collaboration measured inclusive cross sections of 7r+ and 
n~ meson production in the interactions of 31 GeV/c protons with carbon 
nuclei at forward emission angles (0 - 420 mrad). The collaboration also pre- 
sented predictions of Monte Carlo models - FLUKA, VENUS and UrQMD, 
in comparison with the data. A careful analysis shows that deviations of the 
FLUKA and VENUS predictions from the data have different tendencies. The 
worst description of the data was observed for the UrQMD model results. 

All the models assume the creation of quark-gluon strings in the interac- 
tions, but it is complicated to analyze the models in order to find the source of 
the deviations. Thus, the quark-gluon string model - FRITIOF (FTF) - was 
implemented in the GEANT4 toolkit and is used to understand the deviations 
mentioned above. It was found that the most important factor influencing the 
FTF calculations is the sampling of quark-gluon string masses. The other fac- 
tors/parameters are not essential for a description of the data. Also, a good 
description of the data is achieved by the FTF model. 

Some selected results of the NA61/SIIINE collaboration 1 in comparison with the model predictions 
obtained by the collaboration, and taken from Fig. 23 of Ref. [1 together with its figure caption, are 
reproduced in Fig. 1. 
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Figure 1: Laboratory momentum distributions of tt"'" mesons produced in production p+C interactions at 
31 GeV/c in different intervals of polar angle (9). The spectra are normalized to the mean 7r+ multiplicity 
in all production p+C interactions. Error bars indicate statistical and systematic uncertainties added in 
quadrature. The overall uncertainty (2.3%) due to the normalization procedure is not shown. Predictions 
of hadron production models, Fluka2008 (solid line), Urqmd1.3.1 (dashed line) and Venus4.12 (dotted 
line) are also indicated. 
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It is difficult to estimate the size of tlie model prediction deviations from the data given the choice of 
scale and normalization. Though, it can be seen that at small angles the VENUS model ^ overestimates 
the data. At large angles, the model predictions underestimate the data. So, the meson yield as a function 
of angle decreases faster in the model than in the experimental data. 

The FLUKA model [31 underestimates the data at small angles a little bit. At the same time, the 
model overestimates the data at large angles. So, the meson yield versus angle decreases more slowly in 
the model than in the experimental data. 

The UrQMD model 0] predictions are not considered in the present work. A special paper [S| is 
devoted to this subject. 

The FTF model predictions with absolute normalization are shown in Fig. 2. As seen, the model is in 
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Figure 2: Inclusive cross sections of tt"*" meson production. Points are experimental data [T|. Lines are 
the FTF model calculations. 

an agreement with the data at small angles. At large angles, the model underestimates the data. Thus, 
FTF has the same tendency as VENUS. A subject of the present paper is an improvement of the FTF 
model. 




Figure 3: Double diffraction and single diffraction considered in the Fritiof model. 

The FTF model is based on the well-known FRITIOF model [B] (see also Ref. 0). The Fritiof 
model assumes that all hadron-hadron interactions are binary reactions, /ii + /i2 h'l + h'2, where h'^ 
and h'2 are excited states of the hadrons with continuous mass spectra (see Fig. 3). If one of the post- 
interaction hadrons remains in the ground state (h\ + h2— > hi + h2') the reaction is referred to as 
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" single diffraction dissociation" . If neither hadron remains in tlie ground state the reaction is referred 
to as "double diffraction dissociation". The excited hadrons are considered as QCD-strings, and the 
corresponding LUND-string fragmentation model is applied for a simulation of their decays. 

The key ingredient of the Fritiof model is the sampling of the string masses. In general, the set 
of possible final states of the interaction can be represented by Fig. 4, where samples of possible string 
masses are shown. There is a point corresponding to elastic scattering, a group of points which represents 
final states of binary hadron-hadron interactions {hi + /12 ^ ^3 + ^4), lines corresponding to the single 
diffractive interactions, and various intermediate regions. The region populated with the brown points 
is responsible for the double diffractive interactions. In principle the mass sampling thresholds may be 
below the initial hadron masses. The string masses are sampled in the triangular region restricted by the 
diagonal line corresponding to the kinematical limit Mi + M2 = Ecms, and the threshold lines. If a point 
is below string mass threshold, it is shifted to a nearest diffraction line. 




Figure 4: Diagram of final states of nucleon-nucleon interactions. The axes are in GeV. 



In the Fritiof model code implemented in Geant4 we take into account binary reaction final states 
differently from elastic scattering. Weights of the hN elastic scattering and the hN diffractive interactions 
have been introduced which give probabilities of the corresponding reactions. In the case of single 
diffraction we sample the points of the diagrams only on the diffraction lines. 

The probability of the single diffraction (Wsd) at high energies is given as: 

Wsd = D exp{-d Yiab), D = 0.805, d = 0.35, 

where Yiab is the projectile rapidity in the target rest frame. The squared mass of the string created in 
the reaction is equal to 
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where P is the light-cone momentum sampled according to the distribution dP/P. P = E — P^ for the 
hadron originating from the projectile, and P = E + P^ for the hadron originating from the target, in 
the CMS. Too is the mass of the hadron left in the ground state. P^ is the square of the transferred 
transverse momentum, < P^ >=0.15 {GeV/c)'^. 

The probability of double diffraction, Wdd = 1 — '^Wsd, for nucleon-nucleon interactions. In double 
diffraction P~ , the momentum of the hadron originating from the projectile and P"*", the momentum of the 
hadron originating from the target, are sampled independently according to the distribution dW oc dP/P. 
The string mass sampling threshold, mth = toq + 280 Me V. 

The Fritiof model assumes that in the course of a hadron-nucleus interaction the string originating 
from the projectile can interact with various intra-nuclear nucleons and thus reach highly excited states. 
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The probability of multiple interactions is calculated in the simplest approximation. The cascading of 
secondary particles is neglected as a rule. Due to these assumptions, the original Fritiof model fails to 
describe nuclear destruction and slow particle spectra. In order to overcome these difhculties we enlarge 
the model by including the reggeon theory- inspired model of nuclear desctruction [51 12] ■ Momenta of the 
nucleons ejected from a nucleus during the "reggeon cascade" are sampled according to a "Fermi motion" 
algorithm presented in Re f. 10 . 

The main parameters of the model are the transferred transverse momentum and the probability of 
single diffraction at a given energy. The form of the light-cone momentum distribution can be changed 
also. 

It was checked that decreasing (increasing) < > leads to an increase (decrease) of the inclusive 
cross section maxima at 6 < 140 mrad. There is no effect at larger angles. 

Variations of Wsd change the integrals of the cross sections without changing their forms. 

The last possibility - a change of the form of the light-cone momentum distribution, turns out to be 
successful. Calculations performed assuming a flat P distribution {dW oc dP) for the double diffraction 
dissociation processes are shown in Fig. 5 by solid lines. As seen, they are rather close to the experimental 
data. 
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Figure 5: Inclusive cross sections of tt"*" meson production. Points are experimental data Pl. Solid and 
dashed lines are the FTF model calculations with flat and 1/P distributions, respectively. 



It is obvious that average masses of strings are larger with the flat distribution than with the 1/P 
distribution. Thus, multiplicities of produced particles in the target fragmentation region, at large angles 
and small momenta, are larger in the first case. At the same time, the multiplicities in the projectile 
fragmentation region, at small angles and large momenta, become lower due to kinematics with the flat 
distribution. 

From this point of view, we can suppose that the masses of the strings are larger than needed for a 
description of the data in the FLUKA model. Results of the VENUS model point to low string masses 
in the model. All of the features of the model can be improved. 

Note that until now there were only restricted possibilities to check string mass distributions. In 
this sense, the NA61/SIIINE data are unique. It would be well to continue the experimental studies, 
reach higher statistics, and extend the measurements on K-mesons and protons. The proton spectra are 
extremely interesting because Monte Carlo models give very different predictions for them. 
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Conclusion 

1. To reach agreement of the Geant4 FTF model with the experimental data of the NA61/SHINE 

collaboration, it is necessary to replace the 1/P distribution of the squared string masses by the 
flat distribution in the double diffraction dissociation processes. 

2. The string masses sampled in the FLUKA model are probably a little bit higher than necessary for 
a correct description of the data. 

3. The string masses sampled in the VENUS model are too small. 

The author is thankful to D.H. Wright, G. Folger and V. Ivanchenko for interest in this work and 
important remarks. 

References 

[1] NA61/SHINE Collaboration (N. Abgrall et al.)// Arxiv::hep-ex 1102.0983. 

[2] K. Werner, Nucl. Phys. A525, 501c (1991); 
K. Werner, Phys. Rep. 232, 87 (1993). 

[3] A. Fasso et al., CERN-2005-10 (2005); G. Battistoni et al, AIP Conf. Proc. 896, 31, (2007). 

[4] S.A. Bass et al.// Prog. Part. Nucl. Phys. 41, 225 (1998); M. Bleicher et al.// J. Phys. G 25, 1859 
(1999). 

[5] UrQMD patch 

[6] B. Andersson et al., Nucl. Phys. B281, 289 (1987) 

[7] B. Nilsson-Almquist and E. Stenlund, Comp. Phys. Comm. 43, 387 (1987). 

[8] Kh. Abdel-Waged and V.V. Uzhinsky// Phys. Atom. Nucl. 60, 828-840 (1997); Yad. Fiz. 60, 
925-937 (1997). 

[9] Kh. Abdel-Waged and V.V. Uzhinsky// J. Phys. G24, 1723-1733 (1997). 
[10] 5) EMU-01 Collaboration (M.I. Adamovich et al.)// Zeit. Fur Phys. A358, 337-351 (1997). 



5 



